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TECHNICAL NOTE 4281 


SECOND-ORDER SEENDER-BODY THEORY - 
AXISIMMETRIC FLOW 
By Milton D. Van Dyke 

SUMMARY 


Slender-body theory for subsonic and supersonic flow past bodies of 
revolution is extended to a second approximation. Methods are developed 
for handling the difficulties that arise at round ends. Comparison is 
made with experiment and with other theories for several simple shapes. 


INTRODUCTION 


Slender-body theory is the useful approximation introduced into 
fluid mechanics by Munk (ref. 1) for calculating the lift of airships, 
and extended to slender lifting wings in compressible flow by Jones 
(ref. 2). For 6uch problems concerned with lift, its simplicity Is such 
that the solution is independent of Mach number, and Is found merely by 
solving Laplace's equation in two dimensions. 

The theory becomes only slightly more complicated when the thickness 
of a body is of concern. Then the solution includes a logarithmic term 
proportional to cross-sectional area that varies with Mach number, as 
was shown by Ward (ref. 3) in the case of supersonic flow past general 
slender shapes. The analogous result for subsonic flow was found inde- 
pendently by Keune (ref. 4), Heaslet and Lomax (ref. 5) , and Adams and 
Sears (ref. 6). 

Because slender-body theory is so simple and useful, it is natural 
to attempt to improve Its accuracy by Including nonlinear effects in 
higher approximations. Thus, for bodies of revolution in supersonic 
flow, Lighthlll (ref. 7) found the second-order slender-body solution 
for the crossflow due to incidence, and Broderick (ref. 8) attacked the 
flow at zero angle of attack. Recently Light hill has outlined the second 
approximation for supersonic flow past general shapes (ref. 9) • The only 
application to noncircular shapes is the solution for the elliptic cone 
at zero incidence (ref. 10). These four references constitute the liter- 
ature on this subject, aside from papers by Adams and Sears (ref. 6), 
Legras (ref. 11), and Keune (ref. 12), who ignore nonlinear effects and 
seek only a closer approach to the full linearized solution. 
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The present paper is devoted to second-order slender-body theory in 
subsonic as well as supersonic flow, and is restricted to bodies of revo- 
lution. These are the simplest and most practical shapes, and serve to 
illustrate the methods that will be required later in treating bodies of 
general cross section. Only zero angle of attack is considered because 
Lighthill's treatment of the crossflow at supersonic speeds is entirely 
satisfactory, and could readily be extended to subsonic speeds. On the 
other hand, Broderick's solution for the present problem of zero incidence 
at supersonic speeds is so enormously more complicated than necessary that 
it could probably never be applied to any shape except the cone. 

The formal theory set forth here is relatively simple, being comprised 
in equations (l) to (13 )• Complications appear, however, in the case of 
stagnation points, to which a considerable portion of the paper is devoted. 
It is shown that real difficulties arise only for round noses, and that 
for subsonic flow they can be overcome by comparison with the known solu- 
tion for a paraboloid. Only the region spanned by the body is considered, 
though the flow upstream and downstream could be treated in the same way. 
The second approximation, like the first, depends upon an Integral that 
Is the counterpart for slender bodies of revolution of the "airfoil 
integral" of subsonic thin-wing theory (ref. 13 ). 

This investigation was begun in 1953 > inspired by a suggestion of 
Max. Heaslet, to whom the author is indebted also for subsequent helpful 
discussions. Some of the main results were presented at colloquia at 
the University of Manchester and Fort Halstead In 1954 and 1955* Comple- 
tion has been delayed by the search for a method of treating round noses, 
which was only recently found (ref. 14). 


FORMAL SECOND APPROXIMATION 
Resume of Second-Order Problem 


Consider a uniform subsonic or supersonic stream flowing paBt a 
slender body of revolution at zero angle of attack (sketch (a)). The 


r 



question of just how smooth and 
slender it must be will be con- 
sidered later, but the nose (and, 
in subsonic flow, also the tail). 
If not pointed, is assumed to be 
x no blunter than round. 


Vorticity affects the flow ‘ 
Sketch (a).- Notation for body of only in the sixth approximation, 
revolution. and below that the velocity dis- 

turbances induced by the body 
(referred to the speed U of the free stream) are the gradient of a 
perturbation potential $. Linearized theory is concerned -with a first 
approximation cp that satisfies the Prandtl-Glauert equation 
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(1 -M 2 )cp xx + cp rr + -^ = 0 (1) 

(Principal symbols are defined in appendix A. ) If one attempts to improve 
the linearized solution, the second approximation 0 must satisfy the 
iteration equation (ref. 15) 


(l-M 2 )0 xx + 0 rr +%= M 2 


y +1 


M 2 + l 


- M 2 ^ x 9xx ' 


+ 2cp r cp xr + cp r . 2 cp 


rr 


( 2 ) 


The boundary conditions are that the perturbation potential vanish 
radially far from the body (actually at the bow wave in supersonic flow), 
and that the flow be tangent at the surface. To first- and second-order 
accuracy, this tangency condition is 


« 

II 

8- 

at r = R(x) 

(3a) 

0 r = (l+9x)R’ 

at r = R(x) 

(3b) 


With the velocity potential determined, the pressure coefficient is 
given to second order by 


C P = “20X - 0 r 2 - (1 - M 2 )cpx 2 + ^cp^r 2 + £ M 2 ^ 4 (^ ) 

In the slender-body approximation the first term in equation (l) 
can be neglected, except insofar as It appears in the distant boundaiy 
condition. Similarly, for second-order slender-body theory, various 
terms in equation (2) can be omitted (ref. 9)- However, this simplifica- 
tion is unnecessary here because a particular integral of equation (2) 
itself is known; and it would actually complicate the distant boundary 
condition. 


Resume of First-Order Slender-Body Solution 


Slender-body theory is a further simplification beyond linearization 
that describes the flow only in the immediate vicinity of the body - more 
precisely, within a distance from the axis of the order of the local body 
radius. It can therefore be extracted from the linearized solution by a 
limiting process. Similarly, the second-order slender-body theory sought 
here represents the first two terms of an asymptotic series, and can be 
extracted from the full second-order solution. 

For the first-order slender-body solution we adopt the procedure of 
Keune (ref. 4) and Heaslet and Lomax (ref. 5) as being simpler and more 
physically appealing than the methods of Fourier and Laplace transforma- 
tion. The appropriate solution of the linearized equation (l) that 
vanishes far from the body is 
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r F(e)dg 

n /( x -|) 2 + g 2 r 2 


3 = \/i - m 2 for subsonic flow 


cp(x,r) = 


\x-Br 


F(|)d| r - g— — 

* — " , B = vMm -1 for supersonic flow 

- B 2 ^ 


This may be regarded as resulting from a distribution along the axis of 
the body of sources of strength proportional to a function F(x) that is 
to be determined from the tangency condition. Differentiating and inte- 
grating with respect to x gives 1 


cp(x,r) = 




J F(i)cosh _1 ~| 


Then approximating asymptotically for small r in the integrand (and 
also in the upper limit for supersonic flow) gives , near the body. 


~ t <lc S F ^ sgn ( x “^ 2n = ]? ( x ) Zn F(£)sgn(x-S)ln|x-||d£ 

a S (7a) 

"I?/ F(l)Zn = F ( x )* n P(l)ln(x-6)dE (Tb) 


This is the result of Heaslet and Lomax (ref. 5). 

Alternative forms of the integrals that are a great deal simpler for 
either analytical or numerical evaluation were given by Schultz-Piszachich 
(ref. 1 6 ). Excluding an infinitesimal neighborhood of the point x = | 
from the range of integration, carrying out the differentiation indicated 
in equations ( 7 ), adding and subtracting a logarithmic term, and then 
allowing the excluded neighborhood to vanish leads to 


1 If the body has pointed ends (so that P = 0 there), the procedure 

cun be simplified, and it is only necessary to integrate by parts. How- 
ever, we contemplate treating round ends also, at least in subsonic flow. 
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r 


F(x)2n 


pr 


*/ 


F(x) -P(l) 


<p(x,r) =< 


2\/(x-a) (b-x) 

Br r X F(x) -F(£) 

F(X)1 “ 2^ + i x - t - ' dS 


x-i 




(8a) 

(8b) 


The superiority of* these forms is clear if F(x) is a polynomial; then 
the integrands in equations (8) are simply polynomials , -whereas those in 
equations (7) contain logarithms. 


Imposing the first-order tangency condition of equation (3a) new 
determines the source strength F(x) in terms of the body radius K(x) as 

F(x) = R(x)R« (x) 


Hence the first-order slender-body solution has the form 

cp = F(x)2n r + G(x) 


where 


F(x) = R(x)R> (x) 


F(x)2n 


G(x) = < 


2,J(x-a) (b 




F(x) -F(g) 


x-i 


d£ 




(9a) 

(9b) 

(9c) 

(9a) 


The pressure coefficient on the surface of the body is given by 

C-n = -2[F’ (x)2n R(x) + G' (x)] -R ,2 (x) (10 ) 

"e 


Second Approximation 


The slender-body solution of equation (9a) Is clearly a solution of 
Laplace's equation in the cross plane, -which is the Prandtl-Glauert equa- 
tion (l) with the term (l “M 2 )^^ omitted (except insofar as it is 
implicit in the boundary condition far from the body). This linear term 
must therefore be taken into account in the second approximation in addi- 
tion to the nonlinear terms on the right-hand side of the iteration 
equation (2). Hence the slender-body Iteration equation is 
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0 r r + T' = “ C 1 " M 2 }<Pxx + ^ 2 (^T~ ^ + 1 " M^^xx + 2c PAr + ( P r 2( Prr 

(11) 

A particular integral for the linear term on the right, which 
vanishes far from the holy, is given by 2 

- ^(1 -l)r 2 (F"ln r + G" -F u ) 

A particular integral for the nonlinear terms is known to be given by 
(ref. 15 ) 

„ 2 r r 7+1 M 2 \ i 3 " 

- M 2 (F'ln r + S')(p In r+G* 2 -^ WZ\ |§J 

The complete second approximation is the sum of these two particular 
integrals, plus a complementary solution that will have just the form 
of the first approximation, equations (9). Hence the second-order 
slender-body solution for the perturbation potential is 

0(x,r) = (F + f)ln r + (G + g)+^(lf -l)r 2 (P"ln r + G"-F”) + 

M 2 (F'ln r+G')(F In r + G + ffi 1 )-^ (12a) 

where 

N = -n = 2L±i - ~~ (12b) 

Here f(x) is the second-order increment in source strength. Imposing 
the tangency condition of equation (3b) determines it as 

f(x) = (l - 2M 2 )FF’ In R - M 2 MFF ! + (l - M 2 )FG-’ - l^F 1 G - 

| M 2 || + | (1 -M 2 )R 2 (F"lnR + G" F") (lSc) 

2 This result can also be obtained by retaining secondary as well as 
leading terms of the expansions in any of the conventional derivations 
of slender-body theory. If the Heaslet-Loinax method is followed, it is 
necessary to differentiate and integrate with respect to x, as in going 
from equations (5) to (6), two more times in order to avoid divergent 
integrals . 
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K 


Then g(x) is related to f (x) in the same -way that G(x) is to F(x): 

r 


f (x)2n 


S(x) =< 


2 \l (x-a)(h-x) 


if 


! x 


d| 


f(x)ln b + f*tL * L z J & ± u 

2(x-a) J x-| 


(124) 

(I2e) 


However, it will be seen later that this formal result fails at round 
ends. The proper expression for g for round ends will be given in 
equation (40) for incompressible flow and equations ( 52 ) for subsonic 
flow. 


On the surface of the body the expression of equation (4) for the 
pressure coefficient can be simplif ied, using the tangency conditions of 
equations (3)* to 

Cp s = -20 Xg -R« 2 [l + (2-^)cPx s ] + (^-1 )c Pxs 2 +|m 2 R' 4 (13) 


Slender-Body Integrals 


The second-order, like the first-order slender-body solution, is 
seen to require only the evaluation of the "slender-body integrals" 


i-H-/ 


J x |f(x)| = 


F(x) -F(g) 

N-i| 

r x g(x) -F(g) 


Jo X " S 


d| 


d| 


(subsonic) 

(supersonic) 


(l4a) 

(l4b) 


'a 


and their first three derivatives with respect to x, which involve 
integrals of the same form: 


tV _ 11 (x) - F' (j) F(x) - F(a) F(b) - F(x) 

a J a !x _ |( x-a b-x 

J*' = r X F»(x) -F'(g) { F(x) -F(a) 

a Jo x - | x-a 


( 15 a) 

( 15 b) 


* 


etc. 
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Note that only a single integral Is actually Involved, since the subsonic 
one is related to the supersonic one by 

I*|f(x)| = J x |f(x)| + J b{ F ( x )} > a < x < b (l6) 

However, it is convenient to list both. 

As vith the analogous airfoil integral of subsonic thin --wing theory, 
these integrals can be evaluated analytically for a wide variety of func- 
tions. A short table is given in appendix B. 


Transonic Approximation 


It will be seen in later examples that the analytic form of the 
second approximation is rather complicated even for simple shapes. A 
further approximation that yields considerable analytic simplification 
and a remarkably elegant result is that of the transonic small-disturbance 
theory. In that approximation one retains of the nonlinear terms only 
the one that is dominant near Mach number unity, so that the full equa- 
tion of motion is simply 

d-M 2 )0 xx+ 0 rr+ ^= (r + D^ (it) 

In plane flow the accuracy is improved by keeping a factor M 2 in the 
right-hand term, but a test with the exact solution for cones in super- 
sonic flew suggests that the advantage does not persist in axisymmetric 
flew. The effect of retaining the M 2 is simply to change (y+l) to 
M 2 ( 7 +l) in all that follows. 

If one attempts to solve this simplified equation by iteration on 
slender-body theory, the second approximation 1 b, from equations (12), 


0(x,r) = (F + f)2n r + (G + g) + 


7 + 1 

2(M 2 -1) 


F(F’ In r + 0* ) 


( 18 a) 


where 


F = HR' 


f = - 


7 + 1 


2 (M 2 - 1 ) 


FF* 


( 18 b) 

(18c) 


and G(x) and g(x) are related to F(x) and f(x) by equations (9c) and 
(9d) and (l2d) and (l2e). 
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This result can he simplified because the two second-order terms 
in In r cancel. The second-order increment in velocity potential is 
thus found to be a function only of x, given by 



M <1 
( 19 a) 

( 19 b) 


Here a correction for round ends that will be derived later (eqs. (52)) 
has been included in the subsonic case as the last two terms in the 
bracket of equation (l9a). A corresponding correction should probably 
be found also for the supersonic case; if so, equation ( 19 b) does not 
apply to a round nose and the last term might as well be omitted. 


This incremental potential may be regarded as representing a plane 
wave st andi ng normal to the body axis whose amplitude is independent of 
radius within the slender-body approximation (although it of course 
attenuates at distances large compared with the local body radius, where 
that approximation fails). It can be shown that this result holds also 
for bodies of general cross section, where F(x) in equations (19) is 
A' (x)/2jt if A(x) is the cross-sectional area. 

The surface pressure coefficient is given, in the approximation of 
transonic small-disturbance theory, by 

C pg = -20 x - R ' 2 (20) 

Near round ends in subsonic flow the first-order pressure coefficient is 
infinite like x*" 1 , and the second-order increment like x~ s (so that 
neither is integrable for drag). The same is true for a round nose in 
supersonic flow (except that, as just remarked, the second-order increment 
given above may not be even formally correct ) . 


If 

and the 


round ends are excluded, the drag in supersonic flow is integrable, 
second-order increment is given by 


_^gL_ = 2 * -2 1±1 

M 2 -l 


r-cb) 




1 

2 


rV b F(x) -F(g) 

4 4 


F r (x)F' (| )df dx 


J 


(21a) 


Just as in plane flow past an airfoil, this may be either positive or 
negative, according as the body is fatter near its nose or tail. For 
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if round tails are also excluded, setting F(x) -*• F(-x) shows that 
reversing the direction of flow changes the Sign of the drag increment 
of equation (21a) . This means , in particular, that the supersonic drag 
increment is zero for a body with fore-and-aft symmetry. 

The corresponding expression for subsonic flow (with round ends 
excluded so that the drag is integrable) is 


1 

2 


AgD 


P 

'co 


Jt 


7 + 1 
1 -M 2 


b nb 



-'a 


F(x) -F(£) 

( x “ M 


F'(x)F'(|)d£ dx 


( 21 b) 


This differs from the second term in equation (21a) only in having the 
absolute value signs. As a consequence, however, it can be shown that 
this drag increment is zero in conformity with D'Alembert's principle. 


Osvatitsch and Bemdt have shown (ref. 17) that the transonic small- 
disturbance approximation together with the slender-body approximation 
implies a similarity rule for surface pressures on affinely related 
axisymmetric bodies of thickness ratio t, according to which 



+ 2 n(T 2 jl -M 2 |) 


= P 


M 2 -1 ~ 
(7 +l)t 2 _ 


( 22 ) 


Here P is some function of the transonic similarity parameter 
(M 2 -1)/( 7 +1)t 2 . The present theory gives the first two terms in an 
asymptotic expansion of the function P for large values of its 
argument . 


EXAMPLES IN SUPERSONIC FLOW 
Restrictions on Body Shape 


It is implied in the slender-body approximation (as in linearized 
theory) that the velocity disturbances induced by the body are everywhere 
small. This imposes serious limitations on the smoothness of the body, 
even in the first approximation. Thus for bodies of revolution not only 
must the slope R' be small and continuous, but the curvature R" as 
well. Supersonic noses must be at least pointed (R [ sma ll ), and 
supersonic tails and subsonic noses and tails must actually be cusped 
(R' = 0 ). 3 

It Is well known in thin-wing theory that the conditions for 
linearization may be violated locally without destroying the validity 

3 These restrictions are somewhat more severe than those suggested 
by Ward (ref. 2l). He permits discontinuities in curvature and pointed 
subsonic ends, but it is readily verified that these both lead to 
logarithmically infinite surface pressures. 
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of the solution as a whole. This is true also of the first-order slender- 
body approximation. One can permit discontinuities in curvature and even 
slope, and pointed and even round suh sonic ends, provided one attaches no 
significance to the result close to the resulting singularity in pressure, 
or subsequently corrects the solution locally by techniques that have been 
developed for discontinuities in slope in supersonic (ref. l8) and sub- 
sonic flow (ref. 19) and for subsonic ends (ref. 20 ). (Round supersonic 
noses can probably also be permitted, and could be corrected locally if 
the exact solution were known for supersonic flow past a paraboloid of 
revolution. ) 

In the second approximation of slender-body theory (just as in thin- 
wing theory) the restrictions become more severe, and the remedies cor- 
respondingly more complicated, and it Is no longer always true that the 
formal solution is correct except locally. These difficulties are 
greater for subsonic flow because not only are the bodies of interest 
usually blunter (round noses being the rule), but also the restrictions 
are greater (pointed noses being excluded, whereas they are admitted in 
supersonic flow since no stagnation points appear). 


Consequently, application of the present theory to examples of 
subsonic flow will be postponed until nose corrections have been dis- 
cussed. To illustrate the theory, a few examples will now be given for 
supersonic flow. No difficulties appear if the ends are pointed, the 
meridian curve is elsewhere analytic, and one does not inquire too 
closely into the details of the flow near a pointed tail - where the 
flow is actually subsonic and, in any case, the real flow is determined 
by viscosity. 


Cone 


Consider a cone whose slope 
is 5 (sketch (b)), so that the 
body is described by R = 8x. 

With the origin of coordinates at 
the vertex (a = 0), the slender- 
body potential of equations ( 9 ) 
is 



\ 

Sketch (b).- Supersonic flow past 
cone. 


cp(x,r) = 5 2 x(^Zn §J + l) 


(23a) 


Then equation (12c) gives 


(2M 2 ~l)Zn ^+M 2 N + M 2 +| 


f(x) = -S 4 x 


(23b) 
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and equations (12) give as the second-order perturbation potential 
0(x,r) = S 2 x(ln g + l) + S^M^n 2 g-j^(2M 2 -l)Zn ^ + g+l)- 

M^CN+lJ-iB 2 ^jfj- (23o) 


The surface pressure coefficient is, from equation (13 ), 


C * B - K 2 ln ^ ' 1 ) +b4 [ 3b2 “ e sk 



Sketch (c).- Pressure on cone of 5° 
semivertex angle (7 = 1.405). 



M 


Sketch (d).- Pressure on cone of 15° 
semivertex single (7 = 1.405). 


(5M 2 - l)ln Jjr + 2M 2 N MF+IJ 

(23d) 

which is the result first given by- 
Broderick (ref. 8). 

Broderick has compared the 
first- and second-order slender- 
body solutions with the exact 
results (ref. 22) for various cone 
angles. Two cases are reproduced 
in sketches (c) and (d), and the 
second approximation is seen to 
yield considerable improvement for 
moderate cone angles at speeds 
below the hypersonic range. Also 
shown are the results of the 
second-order theory which does not 
involve the further approximation 
of slender-body theory (ref. 15). 

The slender-body simplification is 
seen to reduce the numerical accu- 
racy at high Mach numbers. The 
reason is that, roughly speaking, 
linearized theory and its second- 
order counterpart assume only that 
the thickness ratio t is small, 
whereas the slender-body approxima- 
tion implies also that Bt is 
small- (ref . 23). The latter is a 
more serious restriction at Mach 
numbers appreciably in excess of 
n/ 2. In the subsonic range, on the 
other hand, (3 cannot exceed 1, so 
that the slender-body simplification 
does not significantly reduce the 
numerical accuracy. 
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The corresponding result in the transonic small-disturbance approxi- 
mation is found from equations ( 9 ) and (l 9 }, or simply by discarding al 1 
second-order terms in equations ( 23 ) except those involving N and set- 
ting M = 1. The result for surface pressure coefficient is, in the 
similarity form of equation (22), 


— 1 + 2 Zn(B8) 
S 2 




(24a) 


This series has been extended to a third approximation in unpublished 
work, giving 


5 


—| + 2 Zn(BS) 


= (2 Zn 2 - 1 ) + 



r(7+i)5 2 i 

\12 4/ 

L B 2 j 




(24b) 


Exact numerical solutions of the transonic small -disturbance problem have 
been calculated by Oswatifcsch and Sjodin (ref. 24). The comparison of 
these results shown in sketch (e) gives an idea of the extent to which 
the present theory can penetrate 
into the transonic range. 

As indicated in sketch (e), 
detachment of the bow shock wave 
and attainment of sonic flow just 
behind the shock are both associ- 
ated with a specific value of the 
transonic similarity parameter. 

However, this is not true (in con- 
trast to plane flow) of the "upper 
critical Mach number" at which 
sonic flow is attained at the sur- 
face. This means that the limit 
of convergence of the small- 
disturbance series (such as 
eq. (24b)) cannot be associated 
with the first appearance of a 
subsonic zone in an otherwise 
supersonic flow field. 



(7*l)8 a 

Sketch (e).- Correlation by transonic 
similarity rule of pressure on 
cone; 7 = 7 / 5 - 


Parabolic Spindle 


The analytic form of the second approximation grows complicated for 
shapes other than the cone, except in the further approximation of tran- 
sonic s maJ i -disturbance theory. 
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parabolic spindle. 


The analytic complexity of the 
full second approximation will be 
illustrated for the spindle formed 
by revolving a parabolic arc 
(sketch (f)). Let it be of unit 
length , 'with semi vertex angle 6 
(so that its thickness ratio is 
5/2 ) y and choose the origin of 
coordinates at the nose. Then the 
body is described by 


R(x) = 5x(l - x) (25) 

The first-order slender-body solution of equations (9) is 


The second approximation is found using equations (12) together with the 
integrals of appendix B. Rather tedious computation gives as the surface 
pressure coefficient 


& 2 [^2(l-6x+6x a )Zn ^|^y-l+l6x-22x E j + 

(l-2M a )+6(7M B -‘>)x{a-x)-30(5M ? -3)x !! (l-x) E J jl«(x)+| 2n 2 (l-x> jt3B 2 jz-20x(l-x)+72x a (l-x) a Jln ai 
£ (1-5*^ )+2 ( 33^-2^ )x+ { 333-311}^ )x®+2 (271^-333 )x 3 +4 ( 99-7&P )x* + 
|^6(Tl^-4)(l+Ex)-5{5rf , -3)(l-25x-3nc 2 +25x 3 )J(l-x)in(l-x)+a4 a H^i-l5x(l-x)+it7x=(l-x) a J + 

(l? K 5 + 5)+^12(7^-^)ln ^• + 21-9U! E j>: + |l3(l9-32>r s )ln ^ V? - ^jx 2 * 

^220(5^-3)111 ^ + 625-^ H s Jx 3 +|l25(3-5tf ! hn ! ^\nr] x ‘‘} 


BB(l-x) + 


(27a) 


Here L 2 (x) is Euler's dilogarithm, defined by 



*n(l- |) 
I 





(27b) 


Keune and Oswatitsch (ref. 25) have encountered this function in their 
integral equation theory for slender bodies of revolution in transonic 
flow. They give a short table and reference's to further tables of which 
Powell's (ref. 26) is the most useful. In accord with the second-order 
similarity rule (ref. 27), the surface pressure coefficient has the 
general form 


C-n = t 2 P(x;Bt) + t 4 


p x ( ) +m s p 2 ( ) + (7 + 1) p 


P 3< 


(28a) 
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■where r is the thickness ratio. More specifically, it has the form 
appropriate to smooth slender bodies of revolution 



•^(P^n r+P 2 ) +T 4 j^(p 11 + M 2 p 21 )ln 2 T+ (p 12 + W^p 22 )ln t + 

Pis + m2 p 23 + (r + ■ 1 ) p - P 3 ] ( 2Sb ) 


Enormous simplification results from the approximation of transonic 
small-disturbance theory. The second-order effect is then given by 
equation (l 9a) as 


A 2 0 = - 


7 + 1 S 4 / 
2 M 2 - 1 V X 


y x2+y xsj 


(29) 


Hence second-order effects alter the pressure coefficient at any point 
by 


A 2 Cp = -2 S 4 (l " 15x + 62x 2 - 94 x 3 + 47x 4 ) (30) 


which is plotted in sketch (g). Adding the first-order contribution 
gives, on the surface. 




2(l - 6x + 6x 2 )Zn 


B5(l-x) 


- 1 + l6x - 22x 2 


1 + 


J 


7 + 1 
M 2 -1 



-34(2x-1) 2 + 47(2x-1) 


4 


(31) 


Here the second-order term has been rewritten to make clear that it is 
symmetric about the middle of the body, as indicated in sketch (g), and 
so contributes nothing to the drag. 



Sketch (g).- Second-order increment 
in pressure on spindle. 
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Sketch (h) compares this simple result with Drougge’s measurements 
of pressure on a parabolic -arc spindle of thickness ratio l/6 at M = 1.15 
(ref. 28). It is remarkable that the first and. second approximations give 
successively more accurate values in the region of subsonic flow, which is 
of considerable extent because the free-stream Mach number is somewhat 
below the value (l.l8) for detachment of the bow wave. 



I 

Sketch (h).- Pressure on parabolic spindle with t = l/6 at M = 1.15. 
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SUBSONIC FLOW 


It has already "been pointed out that only under rather severe 
restrictions on body smoothness will the second-order slender-body solu- 
tion be ■uniformly valid over the entire surface of the body. In particu- 
lar, it fails at least locally near stagnation points, and these can 
scarcely be avoided in practice for subsonic flows . We therefore con- 
sider to what extent the formal solution breaks down - and how it can be 
corrected - for subsonic flow past a body that has sharp (conical) or 
round, rather than cusped, ends, but elsewhere satisfies the smoothness 
requirements. (Violations of the restrictions elsewhere than at the 
ends - for example, at discontinuities in slope - could be treated by 
analogous methods; see refs. 18 and 19. ) 


Failure at Subsonic Ends 


Just as in plane flow (ref. 20 ) it turns out that the fo rmal, second- 
order solution for a body with stagnation points may have one of three 
degrees of validity: 

1. Valid except near stagnation points where it predicts 
infinite surface speeds 

2. Invalid everywhere, but finite except near stagnation 
points 

3 . Infinite everywhere 

These three cases are successively more serious (and are accordingly 
associated with successively greater blunting), except that the second 
is more insidious than the third because it gives no warning. 

The distribution of these three cases with respect to nose bluntness 
and Mach number is compared in the following table with the corresponding 
results for airfoils. A regular trend is apparent, bodies being at least 
as critical as airfoils, with the one exception of sharp noses in subsonic 
flow. There, however, the difficulties of case 2 do arise in the various 
components of the solution but happen to cancel in the net result. (Fur- 
thermore, the corresponding airfoil problem could be put into case 1 by 
manipulating, by partial integration, the integrals involved in the 
second-order theory. ) 
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Case 

Bodies 
(using 0) 

Airfoils 

Using 0 

Using i|r 

1. Valid except near 
stagnation points 

Sharp, M = 0 
1 Sharp, M > 0 

Sharp, M - 0 
Round, M = 0 

Sharp, M = 0 
Sharp, M > 0 
Round, M = 0 

2 . Invalid every- 
where, but finite 
except at stagna- 
tion points 

Round, M = 0 

Sharp, M > 0 

Round, M > 0 

3. Infinite 
everywhere 

Round, M > 0 

Round, M > 0 



1 Except for this one case , placement in the first category has 
"been definitely established by actual worked examples (using the 
Janzen-Rayleigh method). In this exceptional case, the placement 
is based instead on the absence of algebraic singularities, which 
might have to be modified by the source eigensolutions discussed 
below, from the second-order solution given below for the spindle. 
It would be well, however, to confirm this classification by 
carrying out the Janzen-Rayleigh solution for a conical tip. 

In the first case, local failure occurs because the true speed is 
proportional to z e near a sharp nose and to near a round nose 

(where z is the distance into the nose and e Is proportional to the 
body thickness), but the slender-body expansion forces these into the 
formal series 


z e = l+s 2 Zn z + 0(e 4 ) 




which are not uniformly valid near z = 0. Recognition of this source 
of the singularities permits one to formulate simple rules for rendering 
the formal solution uniformly valid, with the aid of the correct solu- 
tion for some simple body having the same nose shape (refs. 20 and 13). 

In the second case listed, the over-all failure results from singu- 
lar eigensolutions - extraneous solutions that satisfy the second-order 
equation and the slender-body boundary conditions. They enter because 
of the inexactness of the slender-body tangency condition near the nose. 
The eigensolution is a point source located at the stagnation point if 
one works with the velocity potential (and a dipole if one works with 
the stream function). In plane flow there are at least three simple 
ways to exclude false eigensolutions, but unfortunately none of them is 
applicable in axisymmetric flow. First, the source eigensolution can be 
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excluded in plane flow by working with the stream function (which imposes 
a condition on mass flow that would he violated by an extraneous source). 
However, Stokes' stream function, which should exclude source eigensolu- 
tions in the same way for axisymmetric flow, fails for other reasons to 
yield the correct second approximation (ref. 29). Second, the source or 
dipole eigensolution can simply be deleted as inadmissible in plane flow, 
and the remainder is the correct solution. In axisymmetric flow, however, 
the true slender-body solution may contain a term indistinguishable from 
an eigensolution. Third, there exists a similarity rule that relates 
surface quantities on a single plane airfoil in subsonic flow to those in 
the corresponding incompressible problem (ref. 13 ) , which is free of 
eigensolutions . No such rule exists for bodies of revolution, however, 
nor does the difficulty disappear at zero Mach number. Indeed, it is 
only for round noses in incompressible flow that eigensolutions arise 
(see preceding table); and they can therefore be handled by comparison 
with the known solution for incompressible flow past a paraboloid of 
revolution. 

In the third case listed, divergent integrals arise in the second 
approximation. They can be assigned a finite interpretation only by 
solving the problem by another approximation - either the Janzen-Rayleigh 
expansion In powers of M 2 , or the full second-order theory without the 
slender-body approximation. The Janzen -Rayleigh solution is uniformly 
valid near the stagnation point, and the second-order slender-body solu- 
tion can be extracted from it using the second-order similarity rule 
(ref. 27). The full second-order solution involves source eigensolutions, 
but they can be eliminated by requiring conservation of mass within a 
large contour that lies everywhere far from the region of nonuniformity 
at the nose. (This cannot be done with the slender-body solution because 
it Is not valid far from the body. ) Both these procedures have recently 
been carried out for the paraboloid of revolution, with identical results 
(ref. l4). It will be shown here how any other round-nosed body can be 
treated with the aid of that solution. 


Sharp Ends - The Parabolic Spindle 


Sharp -ended bodies in subsonic flow have stagnation points, but the 
formal second approximation, like the first, is correct except very near 
the tips (case 1 of the preceding table). It can be corrected even there 
by simple rules (ref. 20). However, the region is so minute (being of 
exponentially small order in the body thickness) that the correction is 
usually of no practical significance and will be ignored here. 
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Sketch (i).- Subsonic flow past 
parabolic spindle. 

B(x) = 


As an example, consider again 
the parabolic-arc spindle. Because 
symmetrical bodies in subsonic flow 
induce symmetric disturbances, it 
: is convenient to choose the origin 
at the middle (sketch (i)). The 
spindle of thickness r is 
described by 

(1-x 2 ) (32) 


The slender-body solution of equations (9) is 


9 


T^X 


2(1 - x 2 )ln 


2’JT z ~: 

pr 


:2 

3 +T x 


(33) 


The second approximation is found from equations (12) and appendix B. 
The result for the streamwise velocity component on the surface is 


^ 1+1 * (1 - 3 K 2ta ^- 3 ) 


“3 + 


‘{p 2 (5- 




2M a (l-12x s +15x 4 In 2 “f|) + 


M 2 ^) - (271-325M 2 ~ ^ 


15(3-5^ )x 2 " (27-41M 2 ) 
2 


x In 


1+x 

1-x 


In 


PtaJi-x 2 


1049 5629 


x 2 + 


8 


24 


M 2 +4714 2 n )x 4 


(34) 


The pressure coefficient can be calculated from equation (4), and has 
again the form of equations (28). The maximum velocity, which occurs on 
the middle of the spindle, is given by 


u max 

U 




3 ) 


3 + T< 


10j3 a In 2 


_2_ (%l 6l 
Pt " 


M= 


n 


_2_ 


+ 


jt£ 

12 


(5 - 7M 2 ) M 2 + 3M 2 n 


(35) 
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Again, enormous simplification results from the approximation of 
transonic sma ll -disturbance theory. The second-order effect is then 
found from equations ( 19 } to he the same in subsonic as in supersonic 
flow. It is therefore given by equations ( 29 ) and ( 30 ), with x 
replaced by (l+x)/2 because of the difference in coordinates, and e 
by 2 t. Hence the surface pressure coefficient is given by 

Cp = 2 t2^(3x2 - 1)^2 Zn - 3 ) - 2x2 J - t*( 3 - 34x2 + 47x*) 

(36) 

The spindle in subsonic flow has been treated in the transonic 
s mai 1 -disturbance approximation also by Keune and Oswatitsch (ref. 25), 
who solve an approximate integral 
equation numerically. Their rest 
for the perturbation velocity on ; 
l4.6-percent-thick spindle at 
M = 0.90 (which is nearly the 
critical Mach number) is shown in 
sketch ( j ) to compare reasonably 
well with the present result. 4 I] 
particular, their curve crosses 
that of linearized theory twice o: 
each half of the body, as the 
second-order solution does 
(sketch (g)). 

Drougge (ref. 28) has tested 
a parabolic 6pindle truncated by 
a support sting (cf. sketch (h)). 

If the base lies at x = b 
(sketch (i)), the first-order 
slender-body solution gives 

Cp = 2 t 2 ( 3 x 2 -l)zn -g— ^ -g + f (l+b2) - 3 ( 1 -h)x - 11 x 2 +b 

L P T (l+x)(l-x) c D X J 

(37) 

The algebraic singularity at the comer (x = b) should be corrected by 
the techniques of references 19 and 20. In the second approximation the 
comer introduces divergent integrals (just as a round nose does). This 
difficulty has been avoided by using the second-order increment for the 
complete spindle (the second term in eq. ( 36 )), which should be a 

4 Keune and Oswatitsch solve equation (17) with (7 + 1 ) replaced by 
M*(l -yP)/(l -M*), where M* 2 = (7+ l)M^/[2 + (7 - l)M 2 ] ; this change has 
therefore been made also in equation ( 36 ) in calculating the curve in 
sketch ( j ) . 
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satisfactory approximation away from the comer. The result is compared 
in sketch (k) with the measured pressures at M = O.85, and the second- 
order terms are seen to improve the agreement. 



Sketch (k).- Pressure on parabolic spindle with r = l /6 at M = O.85. 



Incompressible Flow Past Paraboloid 

Consider now the case in which 
x eigensolutions may invalidate the 


second approximation everywhere. 
According to the preceding table, 
this case (case 2 ) can occur only 
for round noses in incompressible 
flow. We consider, therefore, first 
the prototype of round-nosed bodies, 
a paraboloid of revolution. With the 
nose at the origin, it may be 
described by y = -y/ 2 px, where p is 
the nose radius (sketch (l ) ) . Although the infinite paraboloid has 
properly no thickness ratio (or is an ellipsoid of zero thickness ratio), 


Sketch ( l ) . - Notation for 
paraboloid of revolution. 
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Jp formally assumes that role. From equations ( 9 ) the first-order 
slender-hody potential is found to "be, aside from an irrelevant constant 
(which includes the "infinite constant" -(p/ 2 )ln b of eq. ( 9 c)) 

<P = \ P *n ^ (38a) 

Then from equations (12) the formal second approximation is found to be 


Mir 


p In 


r 2 1 

x "8 p 


r 2 

X 2 


( 38 b) 


Here the second-order term is actually incorrect. The exact perturbation 
potential for the paraboloid in incompressible flow is known (e.g., from 
separation of variables in parabolic coordinates) to be 


® \ p in 2 


+ r 2 - [ x - g P 


- I 0 !n f + B -S) + °(P 3 ' » r4) 


(39) 


Thus the formal solution of equation ( 38 b) is seen to be in error by a 
term p 2 /4x, which affects the pressures everywhere. 


This term is an eigensolution for the slender-body problem, because 
it satisfies trivially the equation cp ™ + cp r / r = 0 without affecting the 
slender-body tangency condition of equation (3a). Moreover, it has the 
proper behavior at infinity, because it is in fact the slender-body 
representation of a point source located at (or within a distance of 
order p of) the origin. Thus the exact perturbation potential for a 
point source of strength p 2 /4 located on the axis at x = kp is 



J(x - kp )‘ 


_ t 


+ r 2 


4 x +0(p3 > 


P 2 r2) 


Alternatively, the eigensolution may be regarded as representing a 
second-order uncertainty in the location of the nose. For replacing x 
in equation ( 38 b) by x - (p/2) yields the correct result of equation (39). 


Eigensolutions at Round Ends in Incompressible Flow 


The extraneous eigensolution arises in the formal solution for the 
paraboloid because of the inexactness of the tangency condition near the 
nose; consequently, just the same error will arise for any other body 
having a round nose of the same radius. That Is, the formal second-order 
slender-body potential will be too small by an amount 
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Ag = 


Pa 2 

4(x - a) 


where p a is the radius of the nose, located at x = a. A corresponding 
error will arise at a round rear end, where the eigensolution is the 
slender-body representation of a sink rather than a source. Hence, the 
formal solution of equations (12) can be corrected by calculating g(x) 
not from equation (l2d) but from 


g(x) = f(x)ln 


1 

2«J(x -a)(b -x) 








(4o) 


This modification gives a solution that is valid to second order 
except within a distance of the order of the radius from either round 
end, where singularities remain. That is, removal of the spurious 
eigensolution by means of equation (4o) reduces the difficulty from 
case 2 of the preceding table to case 1. For example, the surface speed 
on a paraboloid of revolution is found, either from equation (39) or 
simply from Muhk's rule (ref. 30) that the speed on any ellipsoid sub- 
jected to incompressible flow along an axis is the projection of the 
maximum velocity, to be 


1 

U 



(4la) 


Expanding this formally for small p yields 

V _ -1 _1 P x 3_ Pf _ 
U 4 x 32 x 2 


(4ib) 


and this is also the result of the present theory, the first two terms 
being the usual slender-body result, and the third the second-order 
increment after removal of spurious eigensolutions . The remaining 
singularities are such that even in first-order theory the integral for 
drag calculated from surface pressure is divergent (though this is not 
a serious difficulty because the drag is known to be zero). 


Rules for Rendering Solution Valid Near 
Round Ends in Incompressible Flow 


The singularities remaining at a round nose can be eliminated, and 
a uniformly valid approximation obtained by applying simple rules to the 
formal solution. Derivation of these rules requires a knowledge of the 
exact solution for some body that matches the one under consideration 
near its nose. The paraboloid of revolution is the prototype of round- 
nosed axisymmetric bodies. It was shown in reference 20 that the ratio 
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of the exact solution for the paraboloid (eq. (4la) ) to its formal series 
expansion (eg.. (4lb ) ) serves as a multiplicative correction factor for 
any round-nosed body. This rule renders the solution correct to second 
order 5 for uncambered airfoils (to which it also applies), but only to 
first order for bodies of revolution. 

A second-order rule for bodies was derived (ref. 20) by considering 
the exact solution for an ellipsoid (or hyperboloid) which matches the 
nose more closely than does a paraboloid. It should be pointed out that 
in this case one cannot simply use the ratio of the exact solution to 
its formal expansion because this would introduce a spurious stagnation 
point at the remote end of the ellipsoid. What one actually requires is 
the exact solution for a semi-infinite body that matches the nose to the 
required order, and this can be extracted from the solution for the 
ellipsoid. 

The result is that for a body of revolution having a round nose at 
x — o, described by 


r 2 = R 2 (x) = 2px-Bx 2 +. . . (42a) 

the formal second-order slender-body solution "q " for surface speed 
is converted into a uniformly valid second approximation q 2 by the rule 

0-2 _ 1 f 0-2 f X 3-1 1 . 

u " V U + 2 U "8 ' 

where "q ” is the first-order part of 

A body with two round ends can be treated by applying this rule 
twice, shifting coordinates so that in equations (42) x is always meas- 
ured into the end. The result can be simplified somewhat to the follow- 
ing. For a body having round nose and tail of radii p & and located 
at x = a and x = b: 



r ”<l " 


U 


J(l+A a )(l+?*b) L u 2 

P, 


+ § ( A a +? 'b) 


y 1 


= S 2. TO 

Aa 2(x - a) 4 > 


^b = 


Pb 


2(b-x) 4 


4 B b 


Corresponding rules for treating surface pressure directly have been 
given for airfoils (ref. 13 ), and could readily be deduced also for bodies. 

5 As pointed out in reference 20, the order of terms is counted in 
such a way that disturbances in velocity or pressure are always of first 
order. Thus a first-order term is only of 0(T 2 £n t) near the middle of 
a slender body, but Is 0(1) near a stagnation point. 
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Rules for treating sharp ends have been given in reference 20; but, as 
discussed in connection with the spindle, the region of nonuniformity is 
so small as to be of no practical significance. 


M 


Example: Incompressible Flow Past 

Ellipsoid 

x Consider a slender ellipsoid 

of revolution in incompressible 
flow. With axes chosen as shown in 
Sketch (m).- Ellipsoid of revolu- sketch (m), the ellipsoid of thick- 
tion in subsonic flow. ness ratio t is described by 



r = R(x) = Wl - x 2 

Equations (9) and (12) give 

F(x) = -t2x 

G(x) = -r^x^n 2*^1- x 2 - 1^ 
f(x) = -T^ln f + §)x 


(43) 




(44a) 


The radii of the nose and tail are p = p, = t 2 , so that equation (4o) 
gives 


k(x) = -r 4 ^2n ^ + j^x ^Zn 2^1-x 2 - 1^) r 4 — 


■X2 


(44b) 


Then from equation (12a) the second-order slender-body solution for the 
perturbation potential is 


0 = r 2 1 + ~ +§) x (* n 




T^r^x 


3 -x 2 1 

(l-x 2 ) 2 “2 


T 4 


X 

1 - X 2 


(44c) 


It can be verified that this is the asymptotic expansion, to this order, 
of the known exact solution for flow past an ellipsoid. 

The streamwise velocity component and resultant speed on the surface 
are found to be 
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^ = 1 + r^Zn — ) 

U \ T 1 -X 2 / 


, p 2 1 1 + x 2 2 1 x 2 

Zn 2 Zn r-+- 


T - 1 + t2 


Zn 


T 2 1-3(2 T 4 ( 1 - X 2)2J 

I -I C 1 + r^)J + I- 1 l “ I - b - 


1 1 ,3 1-1 

4 1 - x 2 8 (1 -x 2 ) 2 J 


(45a) 


(^5b) 


Applying the rule of equation (42c), -with P a = = B & =B^ = t 2 , gives 

the uniformly valid result for surface speed (or pressure coefficient) 




C-n — 


MrMM T <§-r^)J' 


(46a) 


The exact result is 


1 + T 2 


£L 

U 


sech -3 ^ - s] 1- t 2 
\l 1- t 2 - T^sech 

(— rST 


(46b) 


As an extreme test, the approximate and exact values are compared in the 
following table for an ellipsoid of thickness ratio t = 1 / 3 : 


±x 

0 

0.4 

0.8 

0.9 

0.95 

O .98 

1 

q/U 

1.122 

1.110 

1.025 

0.924 

0.788 

0.584 

0 

V u 

1.114 

1.103 

1.018 

.918 

cy 

CO 

t- 

m 

.580 

0 


Subsonic Flow Past Paraboloid 


The remaining case to be disposed of is that of subsonic flow past 
a round nose. This is case 3 of the table on page 18, in which the 
formal second approximation leads to divergent integrals. This will be 
illustrated for the paraboloid (sketch (z))j and comparison with the 
known correct solution will again provide appropriate corrections. 
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The first-order slen&er-body solution for the paraboloid is found 
to be independent of Maeh number, so that it is given by equation (38a). 
However, the potential is indeterminate to within an additive constant 
which was dropped there but must now be retained for purposes of the 
comparison. Hence the slender-body potential is written to include an 
arbitrary constant K an 



In the second approximation, equation (l2c) gives as the increment 
in source strength 

f (x) = - | M 2 p 2 J (48a) 

and difficulties appear because this is not integrable at the nose. The 
function g(x) of equation (l2d) may he written formally as 



(48b) 


If e is here regarded as small, all difficulties have been concentrated 
into an integral over a short portion of the nose. The integral diverges 
so that it is meaningless as it stands, nor can any a priori significance 
be assigned to it as a finite part. The proper interpretation is rather 
to be found from comparison with the known solution. 

The formal second approximation of equation (12a) thus becomes 



(48c) 


whereas the correct result has been shown to he (ref. l4), aside from an 
irrelevant constant. 


0 = 


'I* 
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+ H 3 2 f 
r 2 8 


k-r* Vl 

X X 2 / 4 
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+ 2 In 


2x 

Pr' 


»)- 


r 2 J 

(49) 


These two expressions agree if the divergent integral is interpreted 
according to 
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(50) 


Eigensolutions at Round Ends in Subsonic Flow 


Consider now the general case of a body having a round nose of radius 
p a at x - a and a round tail of radius at x = b. The second- 

order increment to source strength f (x) will consist of a regular func- 
tion f*(x) plus the singular terms 


f(x) 



(51) 


These give a divergent integral / d|/| at either end of the body. 

^ o 

Each of these integrals can be interpreted according to equation ( 50 ) in 
terms of the corresponding radius p a or p^ and a constant Kg_ or 
that can be determined from the first-order solution. The result is 
that in place of equation (l2d), the function g(x) is given by 


g(x) = f(x)Zn 


2 n /(x - a) (b -x) 


■ 1 / 


lb f*(x) -f*(|) 


l*-e| 


u + 



») 


+ 2n 


8p a (x - a) 2 
P 2 (t-a) + 



Pb 2 
b -x 



SpfaCb -x) 2 
P 2 (b - a) + 



( 52 a) 


where 




( 52 b) 


K, = lim 
x-*a 


G(x) Zn(x - a) 

— 


K b = lim - — G(x) +J Zn(b -x) 
x-*b[_ Pb tL 

As M tends to zero this reduces to equation (4o). 


( 52 c) 

(52d) 
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As in incompressible flow, this modification renders the solution 
valid to second order except within a distance of the order of the radius 
from either end. (Case 3 of the table on page l8 has been reduced to 
case 1.) The surface speed again contains singularities like (x-a) _1 
in the first-order terms and (x - a) -2 in the second-order terms (cf. 
eq. (4lb) ) . These can be eliminated, and the solution rendered uniformly 
valid, by a rule corresponding to equations (te). Derivation of the 
second-order form of this rule makes use of the formal solution for the 
ellipsoid, which must therefore be found first. 


Example: Subsonic Plow Past Ellipsoid 


Consider subsonic flow past the slender ellipsoid of revolution of 
sketch (m). According to equations ( 9 ), the first-order source strength 
F(x) is unchanged from the incompressible value of equation (Vt-a), and 
G(x) is modified only by insertion of a factor p, so that 


F(x) = -t 2 x 

G(x) = T a x(ln - l) I 


(53a) 


The coefficients of the constant terms in G(x) at the ends are, from 
equations ( 52 c) and ( 52 d), 


Kq, — K*b = 1 + In 


JL 

2^2 


(53b) 


Equation (12c) gives as the formal second-order increment in source 
strength 


f(x) 


T^X 


(2M 2 -l)ln 


_ 2 _ 

pr 


+ M^n - 1) 



+ | !#t 4 


x 

1 -x 2 


(53c) 


The quantity in brackets is the f*(x) required in equation (52a), and 
the remainder is the singular terms of equations ( 51 ) and ( 52 b) that 
lead to divergent integrals. The integral in equation (52a) is trivial 
(being a multiple of the first-order one), with the result that 


g(x) = ~T 4 x 


(2M 2 -l)ln Ip + M^tn -1) -|](ln - 1 ) B 2 ' r 4 5 - 7 *= 


1 4 /, Tn/1-X* 

5 K T \ ln ~T 


+2 1 57^+S ^ 


Zn(l + x) ln(l - x) 
1 +x 1 -x 


(53d) 
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Finally , equation (12a) gives as the second-order slender-hody perturbation 
potential 




1 p 2 -!- 4 — + 1^T 4 X; 

2 1-x 2 


\ Pr l-x 2 A 


2n/1 -x2 


. , A In^^ + n-l)+f t 2 

Pr 1-xVV Pr 


1 

'4 


X 2 

r 2 


( 53 e> 


As a partial check, it may he noted that this result conforms to the 
second-order similarity rule (ref. 27). Differentiating gives the velocity 
components. On the surface of the ellipsoid 


u 

U 



1 

1-x 2 



— + M 2 - l^x 2 
Pt |_ 


3 - x 2 1 1 + x 2 2 
(l -x 2 ) 2 2 1 -x 2 _^ n Pt + 



1 + X 2 
(1 -X 2 ) 2 


In 2 Zn ( 1 + X ) In(l-x) 3 1-^x 2 + 2x 4 ~| 1 x 2 

4(l +x) 2 4(l - x) 2 4 (i-x 2 ) 2 J 4 (i-x 2 ) 2 + 



1 - 5x 2 + 2x 4 l 
(1-x 2 ) 2 J 


(54a) 



1 

2 


1 

1-x 2 




M 2 - M 2 x 2 


3~x 2 

(l“-x 2 ) 2 


1 

2 


1 

1-x 2 


] 




1+X 2 

(1-x 2 ) 2 


• 7 V, o in(l+x) Zn(l-x) 3 l-5xS+2x 4 l 3 11 

4(l+x) 2 “4(l-x) 2 "4 (1-x 2 ) 2 J "8 "4 1-x 2 + 


3 

8 


(1-x 2 ) 2 


+ ^I^n 
2 


l~5x 2 +2x 4 l 

(1-x 2 ) 2 J 


(5^h) 


The maximum speed in the flow field (aside from spurious singularities 
at the ends of the ellipsoid that are to he removed) occurs at the surface 
in the middle, and is 


Imax 

U 





ii-f-i M 2 n 

Pt 2 


1 

4 



( 55 ) 
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ThiB agrees with, the result of Schmieden and Kawalki (ref. 3^-) except for 
the coefficient of M 2 t 4 , which they give as ~5 fb instead of “(^ " 7? 2n 2^. 


They work with the perturbation form of Stokes' stream function rather 
than the velocity potential , which facilitates imposing the condition of 
tangent flow at the body (particularly since they impose the condition 
exactly, and only later extract the slender-body series). However, they 
retain only linear and quadratic terms in the equation of motion. The 
cubic term M 2 cp r 2 cp rr in equation (ll) yields a second-order effect, and 
the same is true of some quartic as well as cubic terms when one works 
with Stokes’ stream function. Furthermore, the linearized equation for 
the stream function is not correct to first order except in the slender- 
body approximation and in any case does not form a proper basis for 
iterating to find the second approximation (ref. 29). Thus for supersonic 
flow past a circular cone, Schmieden and Kawhlki’s procedure was found to 
yield the second-order slender-body solution - correct except for the term 
in M 2 t 4 , This is presumably true in general, so that the disagreement 
in that term found here for the ellipsoid might have been anticipated. 

The present solution predicts a maximum speed slightly higher than 
Schmieden and Kawalki’ s, which does not appear unreasonable in view of 
their comparison with the Janzen -Rayleigh solution to order M s for a 
sphere, in which their speed was somewhat low. 


In the transonic small -disturbance approximation, the surface pres- 
sure coefficient is given simply by 




__1_\ y +1 t 4 1 - 5x2 + 2x 4 
1 -x2/" 2 1 - M 2 (1 -x2) 2 


(56) 


Rules for Rendering Solution Valid Near 
Round Ends in Subsonic Flow 


As in incompressible flow, the ratio of the exact solution for a 
paraboloid to its formal series expansion serves as a first-order 
multiplicative correction factor. Thus it has been shown (ref. lb) 
that the slender-body solution for surface speed is converted into a 
uniformly valid first approximation by the rule 



where Q(2x/p , M) is the speed ratio on a paraboloid of revolution of 
nose radius p at Mach number M. Although Q, has not been found 
exactly, it can in principle be determined to any desired degree of 
accuracy by the Janzen -Rayleigh method. In reference l4 this has been 
done including terms in M 2 , and numerical values have been tabulated. 
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Singularities remain if this first-order rule is applied to the 
second-order solution. A second-order rule is required; and its farm 
can he deduced with the aid of the preceding solution for the ellipsoid. 

Replacing x by x/c in equations (53) to ( 56 ) gives the solution 
for an ellipsoid whose length is 2 c rather than 2 (cf. sketch (m)). 

Then for small values of the distance z = x + c measured from the nose, 
the surface speed is found from equation (54b) to have the form 



1 4- T 2 


C 

+ 


in 


pr 


— 

« 


+ 


T 


4 


M 2 cf. 

TT z 2 


in 


P 2 t 2 c 

2 z 


+ 


(j- + l 

\32 8 









(57a) 


The parameters p and B of equation (42a) are related to the present 
c and -r by p = t 2 c and B = t 2 , and in those terms the above expression 
becomes 





+ 


Mf pf 

14- Z 2 


in 



14- 

0 2 B 




(57b) 


This expansion could be used to form a second-order rule, but the result 
is simplified by first determining the corresponding expansion for a 
general body, and then choosing a simple special case. 


It is clear that, corresponding to equation (57b), the solution for 
any round-nosed body described by equation (42a) will, near its nose, 
have the form 



1 + 




+ 


M 2 pf 
”4" z 2 


in 





( 58 a) 


Here G is a constant, and R x and R 2 are regular functions of z 
between which a relation will now be found. Although equations (57b) 
and ( 58 a) are both singular at z = 0 , their ratio must be regular. 
Dividing the latter by the former and expanding gives 



3^ 


NACA TN 1+281 


1+ R x (z) + (§“§ ln |2b) b + • • .] + {|h x ( 8 )-c + £b §+• • .j- 


and this is regular at z = 0 only if 


0 = nl B+ I Rl(0) 


(58K) 


There is a particular body for which equation ( 58 ) has the simplest 
form, namely that for which, R x (z) = Ra(z) = 0, so that 

"IT = 1 "l& + _T & ln ^ + (^ + i ^ "¥ ^ n ) & B f (59) 

Note that except for the term in B this is Just the second-order slender- 
body solution for the paraboloid. 

Now suppose that the exact solution is known for any semi-infinite 
body having prescribed values of p and B. Then the ratio of it to its 
formal series expansion serves as a second-order multiplicative correc- 
tion factor. Hence a uniformly valid second approximation is given by 

Is exact solution n< l" 2 

U , f P - / il jM 2 P 2 , P 2 P f Z 3 . 2 . rfn\ P £ f 3 „ Ri(Q) ~| P _ , 1 u 

_ 1+ [" ^ + Rl(z) J + lT ^ ^ -2^ + (j2 + f ^ [3 B+ “T-j ^ + R2(z) JJ seml - 

infinite 

body 

(6 Oa) 

Dividing both the series for "q. 2 "/ u and that in the denominator of the 
bracket by equation (59) makes them regular at z = 0, and gives, after 
expansion. 


exact solution 


1+R x (z) +|r 2 (z)+| R x (z)-R x (0) 

M 2 7 ^ 2z 1 3 „e , p 2 3 p P~1 


pl« tl It- Tt 

qg p q.i 

u l+z u 


'semi- 

infinite 

body 


( 60 b) 


This rule may be written finally as e 


S As M-*-Q this reduces not to the rule of equation (k£b) but to an 
equally valid alternative. 
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H'n 11 

, P 

'V + /M 2 1q 2z_l_ 

3 

if +— 

p!,j_ r p] 

Vp ,B ’ J 

u + bz 

u \b 1 p 2 P 32 

8 


z 2 l6 z J 


(6 Oc) 

where Q(2z/p, B; M) stands for the first bracket in equation ( 60 b). 

Note that the terms in the bracket are just those required to cancel all 
singularities in "q 2 "/u. 

A combined rule for two round ends is again found by applying the 
rule twice in succession and simplifying insofar as possible. The result 
is 


Ss. 

u 


2(x - a) 
Pa 



Q 


2(b -x) 
Pb 



1 

¥ 



Pb 

b -x 



P a 2 

4 L(x-a) 2 


2(x - a) 
In — - + 

P 2 P a 


Pb c 


(b-x} £ 


In 


2(b-x) 

P 2 Pb J 


_1_ 

32 





r Pa 2 Pb 2 

3 /^a p a BbPb\t 

_(x-a) 2 (b-x) 2 J 

16 \x - a b - x/J 


( 60 d) 


Mixed Rules Based on the Paraboloid 


The function Q(2z/p , B j M) required In the above rules could in 
principle be dete rmin ed to any desired accuracy by the Janzen-Rayleigh 
approximation . However, the practical details appear almost insurmount- 
able except in the special case of the paraboloid, for which B = 0. 

(The next section shows that even the solution for the paraboloid has not 
yet been carried far enough to yield reasonable accuracy at high Mach num- 
bers.) It is therefore worthwhile to simplify the rules so as to base 
them on the solution for the paraboloid. 

If B is to be replaced by zero in the argument of Q, it must be 
omitted elsewhere in equations (6o). Hence the rule for a single round 
edge, corresponding to equation ( 60 c), is 


U 



M 


ti 

ge 

. U 




’V 

U 




In 

P 2 p 32 8 


M 2 n 


fif i 


(6la) 


Likewise the counterpart of equation (6dd) for two round edges Is 



36 


NACA TN ltafiL 



These rules give a result that is, of course, correct to second 
order except near the ends. It is correct only to first order within a 
distance of the order of the radius p from the ends. Finally, it is 
completely invalid within a much smaller neighborhood of the ends of the 
order of B 2 p (which is proportional to the- sixth power of the thickness 
ratio for a body of unit length). The reason is that the bracket in 
equation (6la) has not been completely freed of singularities, but retains 
a term -3Bp/l6z (which is cancelled in the original rule of eq. (60c)). 
For most practical purposes this distance is so minute as to be altogether 
negligible, and (as indicated by the subscripts) these mixed rules can be 
regarded as yielding a solution valid to first order near the ends and to 
second order elsewhere. 


Example: Uniformly Valid Solution for Ellipsoid 


These rules can now be applied to the formal solution of equa- 
tion (5*t-b) to find a uniform approximation for the ellipsoid. Using 
the combined rule of equation ( 60 d) gives as the uniform second-order 
solution 


q. 2 

f ^ 


^ , i-l * 


Q 


2(1 -x) o If J 2 r 

' 1 > t2; M J{ 1 + T V n ^-2, 


1 (1 -Vp)ln 2 — + 2FPln -^ + M 2 n ^ M 2 j- 


( 62 ) 


The last bracket has of course been freed of singularities, but in addi- 
tion It is seen to be a constant, independent of x. This simple feature 
Is the counterpart, for second-order subsonic flow, of Munk's rule for 
incompressible flow past an ellipsoid, according to which the surface 
velocity is just the projection of the maximum velocity. 
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Using instead the mixed rule of equation (6lb) yields 


U 


= Q 


~2( l + x) 
T 2 


M 


Q 


-2(1 -x) 


s m] |l + ^nl-|y 



(1 - l^)Zn 2 


; + 2 M 2 ln J^ + lVFn- M 2 
Pr Pt 8 2 


3 

8 


1 

1 -x 2 



and the remaining nose singularity is evident. Its effect may he illus- 
trated in the case of incompressible flow. At a distance of one radiu6 
from the nose (which would be 0.02 of the length for a 20-percent-thick 
ellipsoid), equation (63) gives 


U 




The exact solution has - 12 /l 6 instead of -13/16, showing that the result 
is valid only to first (though nearly to second) order. Again, at a 
distance of only Bp = t 4 from the end (0.0008 of the length), equa- 
tion (63) gives 



U 



The exact solution lacks the factor 13 / 16 , but the leading term in the 
pressure coefficient Is nevertheless given correctly as unity, so the 
result is regarded as being correct to first order. This ceases to be 
true only at distances of the order of B 2 p from the end, which is 
0.000032 of the length for a 20-percent-thick ellipsoid. 


Comparison With Experiment 


Matthews (ref. 3 ^) has meas- 
ured pressures over the front half 
of an ellipsoid of revolution of 
thickness ratio l/6 up to Mach num- 
bers of 0.9^0 (the measured criti- 
cal Mach number being 0 . 916 )* The 
ellipsoid was supported from the 
rear by a sting but, according to 
first-order slender-body theory, 
the sting affects the pressure coef 1 
ficient over the front half by less 
than 0 . 003 , which Is negligible. 

The pressures measured at M = 0.900 
are compared with first- and 
second-order theory in sketch (n) . 



Sketch (n).- Pressure on ellipsoid 
of revolution with t = l/6 at 
M = 0 . 900 . 
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Over the middle of the "body the 
with the second approximation, which 
ment over the first. Near the ends, 



Sketch (o).- Speed ratio on parabo- 
loid of revolution at M = 0.90. 


Ames Aeronautical Laboratory 

National Advisory Committee for 
Moffett Field, Calif., May 


experimental values agree closely 
is seen to be a significant improve- 
however, the experimental values lie 
between the predictions of second- 
order theory with and without the 
application of the mixed rule. The 
reason for this is believed to be 
simply that the values of Q, the 
surface speed on a paraboloid, used 
in the rule are inaccurate. They 
were taken from the Janzen-Rayleigh 
approximation including only terms 
in M 2 (table II of ref. l4), which 
is almost certainly inadequate at 
M = 0.900. Indeed, the present 
theoretical results are believed to 
be sufficiently trustworthy that 
one can work backwards to extract 
experimental values of Q, for the 
paraboloid from the measurements on 
ellipsoids. The result is shown in 
sketch (o) in comparison with all 
existing theories. 


Aeronautics 
28, 1958 
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APPENDIX A 
PRINCIPAL SYMBOLS 


abscissa of nose of body 
abscissa of tail of body 

(m 2 - 1) 1/2 

bluntness of nose (see eq. (42a)) 
constant in equations (58) 
pressure coefficient 
half-length, of body 


drag 

source strength in slender-body theory, RR’ 
second-order increment in F 

term independent of r in slender-body potential 
second-order increment in G(x) 
subsonic slender-body integral (see eq. (l4a)) 
supersonic slender-body integral (see eq. (l4b)) 
arbitrary constant 

constant term in slender-body potential (see eqs. (47) and 
(52)) 


length of body 

Euler’s dilogarithm 

free-stream Mach number 

7 +1 M 2 
2 1 -M 2 

7+1 M 2 

supersonic counterpart of n, — — — ^ — — 
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P,Pi,Pa \ 

Pi,P 2 ,P 3 / 

Q 

a. 

R(x) 

Pi^P 2 

r 

U 

u 

X 

z 

p 


general functions in similarity rules (see eqs. (22) and ( 28 )) 

speed on surface of semi-infinite round-nosed "body, referred 
to free-stream speed 

local speed of flew 

radius of meridian curve of "body of revolution 
regular functions 

radius in cylindrical polar coordinates 
free-stream speed 
streamwise velocity component 
streamwise coordinate 

abscissa measured from round end into body 
( 1 - M 2 ) 1 ' 2 


7 

As 

5 

e 

A 

P 

Poo 

T 

$ 

9 

0 


adiabatic exponent 
second-order increment 
initial slope of sharp-nosed body 
small parameter 
(See eq. (42b).) 

nose radius of round-nosed body 
free-stream density 
thickness ratio 
full velocity potential 

first-order perturbation velocity potential 
second-order perturbation velocity potential 
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Subscripts and Superscripts 


( >a 

associated with nose 

c >b 

associated with tail 

( ) B 

value on surface of body 

C K 

first-order value 

( ) 2 

second-order value 

( ) 2;1 

mixed second- and first-order value 

( ^max 

maximum value 

C )* 

regular part (see eq.. ( 52 b)) 

0 

singular part (see eq.. (51)) 

( )’ 

derivative 

C) 

uniformly valid value 

"C )” 

formal value 
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APPENDIX B 


SHORT TABLE OF SLENDER -BODY INTEGRALS 


The integrals appearing in the slender -"body solution were denoted 
in equations (l4) by 


I^F(x) 


hi 


^F(x) -F(g) 


|x-g| 


d| (subsonic) 


(Bla) 


J^jp(x)j- = J' d| (supersonic) (Bib) 

In the supersonic case the notation is designed to emphasize the different 
roles played by x in the integrand and in the upper limit of integration. 
The subsonic integral can be expressed in terms of the supersonic one by 

I^|f(x)| = sgn(x - a) J^-|F(x)j-+ sgn(b - x)J^-|f(x)| (B2) 

of which equation (l6) is a special case. 

For purposes of shifting the origin of abscissas, it is convenient 
to relate the general supersonic integral to that for some standard value 
of the lower limit a, say zero. The desired expression is easily seen 
to be 1 


J?4f(x) 


'} = Jo "^ 1 


F(x + a) 


Combining these last two results gives a useful expression for the general 
subsonic integral in terms of the standard supersonic one; for a < x < b, 

I^F(x)j- = J^“ a |F(x + a)j-+ J^" b |F(x +b)j- (B4) 

A short table of the subsonic and supersonic slender-body Integrals 
is given below. The limits of integration have been taken as a = -1, 
b = 1 for the subsonic case and a = 0 for the supersonic. Results for 


1 The meaning of J^ -a -jF(x +a)j- in conjunction with the following 

table is that one looks up F(x + a) in the column labeled F(x) and in 
the corresponding column labeled J Q replaces x by (x-a). 
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other ranges of integration can he extracted using the two relations 
above (eqs. (B3) and (b 4) ) . Values of the subsonic integral are given 
only for the span of the body, that is, for -1 < x < 1. The table was 
checked where possible using equation (B4). 


F(x) 

H 

«-i 1 
H 


r' a2 ^i aidt 

1 

0 

0 

X 

2x 

X 

X 2 

3x2-1 

2 ^ 

X 3 

y- X 3 -X 

U X 3 
6 x 

X* 

H x *- x =- i 

6 X’ X- 2 

22*4 
12 x 

x^n^l^,. . . 
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2(1-In 2) „/x 
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r 

< 

V l-x 2 (l-In 2s/l-x 2 )-x Bln -l xJ 

— 

In x 

— 

*2 

6 

x In x 

— 

*(*n X+ T _1 ) 

x®In x 



^2 X + f -f) 

In(l-x) 

■B^Oi 

"I In 2 (l-x) 

x In(l-x) 

— 

-xL £ (x) -|xln 2 (l-x)-(l-x)ln(l-x)-x 

x 2 In(X-x) 

— 

-x2L 2 (x) - ^ x?ln 2 (l-x) - 
| (l-x) (l-2x)ln(l-x) - | x - 1 x 2 

x 3 ln(l-x) 

— 

-x3L 2 (x) - | x^I n 2 (l— x ) - 
g (l-x) (2+5X+HX 2 ) ln(l-x ) - 

(M-Hf-) 

x*In(l-x) 

— 

-x*L 2 (x)-|x*In 2 (l-x)- 

(1_x) (f + ^ x + S x2 + S x3 ) ln(1_x) ■ 

(f x+ i x2 + t x3 + f^x*) 

la(l-xs) 

Jt 2 1 _ l+x 

T" 2 ln 1^ 

- - - 

x 2 In(l-x2) 

■ 

— 

x-HnCL-x 2 ) 

(?■ 

'*-* = x= -^§ **) 
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